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Abstract 



This paper is dedicated to the study of the initial value problem for density 
, dependent incompressible viscous fluids in with N > 2. We address the question 

' of well-posedness for large data having critical Besov regularity and we aim at stating 

, well-posedness in functional spaces as close as possible to the ones imposed in the 

(U ' incompressible Navier Stokes system by Cannone, Meyer and Planchon in [fl where 

" ' uq e B/^oo with 1 < p < +00. This improves the analysis of [T3j, [H] and [2] where 

Uq is considered belonging to B^^^ with I < p < 2N. Our result rehes on a new 
a priori estimate for transport equation introduce by Bahouri, Chemin and Danchin 
p in [5] when the velocity u is not considered Lipschitz. 

< 

rH ; 1 Introduction 



> 
(N 
00 



^ ' In this paper, we are concerned with the following model of incompressible viscous fluid 

with variable density: 

' dtp + div{pu) = 0, 
dtipu) + div(pn ^ u) - div{2p{p)Du) + VIl = pf, 
divti = 0, 

. (p,M)/t=0 = (PO,Mo)- 

. Here u = u[t, x) G stands for the velocity field and p = p{t, x) G M"^ is the density, 

I Du = ^(Vn -|-* Vu) is the strain tensor. We denote by p the viscosity coefficients 

' of the fluid, which is assumed to satisfy p > 0. The term VII (namely the gradient 

of the pressure) may be seen as the Lagrange multiplier associated to the constraint 
^ I divTx = 0. We supplement the problem with initial condition {pQ,uo) and an outer force 

^ ' /. Throughout the paper, we assume that the space variable x G or to the periodic 

^ ■ box with period at, in the i-th direction. We restrict ourselves the case N > 2. 

The existence of global weak solution for (II. ip under the assumption that po £ L°° is 
nonnegative, that divuo, and that ^/poUo G has been studied by different authors. It 
is based on the energy equality: 

WVP^mh + \\y^Du{T)\\l2dT = \\^uo\\l2+ J 2{pf ■u){T,x)dTdx. (1.2) 
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Using ()1.2p and the fact that the density is advected by the flow of u so that the 
norms of p are (at least formally) conserved during the evolution, it is then possible to 
use compactness methods to prove the existence of global weak solution. This approach 
has been introduced by J. Leray in 1934 in the homogeneous case (i.e p = 1) and no 
external force. For the non- homogeneous equation (II. ip . we refer to ^ and to [28] for 
an overview of results on weak solution. Some recent improvements have been obtained 
by B. Desjardins in [19], [20] and [21] • 

The question of unique resolvability for (jl.ip has been first addressed by O. Ladyzhen- 
skaya and V. Solonnikov in the late seventies ( see [27j). The authors consider system 
(ll.ip in a bounded domain $7 with homogeneous Dirichlet boundary conditions for u. 

2_2 

Under the assumption that uq ^ W i {q > N) is divergence-free and vanishes on 50 
and that po £ C^(17) is bounded away from zero, the results are the following: 

• Global well-posedness in dimension N = 2 

• Local well-posedness in dimension = 3. If in addition uq is small in W « , 
then global well-posedness holds true. 

Let us mention by passing that for the dimension = 2, O. Ladyzhenskaya and V. 
Solonnikov use a quasi-conservation law for the norm of the velocity and get global 
solutions. 

The case of unbounded domains has been investigate by S. Itoh and A. Tani in [23]. In 
this framework, system (jl.ip has been shown to be locally well-posed. In the present 
paper, we aim at proving similar qualitative results in the whole space or in the torus 

under weaker regularity assumptions. 
Guided in our approach by numerous works dedicated to the incompressible Navier-Stokes 
equation (see z.g [29j): 



{NS) 



dtv + v-Vv- pAv + Vn = 0, 
divv = 0, 



we aim at solving in the case where the data {po,UQ,f) have critical regularity for the 
scaling of the equations and in particular when the initial velocity is in the same Besov 

spaces than Cannone, Meyer and Planchon in [7] for the incompressible Navier-Stokes 

iv_i 

system. It means that we want get strong solutions when uq belongs to Bp^oo with 1 < 
p < +00. By critical, we mean that we want to solve the system functional spaces with 
norm in invariant by the changes of scales which leaves (jl.ip invariant. That approach 
has been initiated by H. Fujita and T. Kato in [23j. In the case of incompressible fluids, 
it is easy to see that the transformations: 

v{t,x) — >lv{ft,lx), V/gM 

have that property. 

For density-dependent incompressible fluids, one can check that the appropriate trans- 
formation are 

{po{x),uo{x)) — > (po{lx),luo{lx)), VZeM. 

{p{t, x),u{t, x),U{t, x)) — > {p{l^t, lx),lu{l^t, lx),l^U{l^t, Ix)). 
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The use of critical functional frameworks led to several new well-posedness results for 
incompressible fluids (see [7], In the case od the density dependent incompressible 

fluids we want cite recent improvements by R. Danchin in [11], [12], [13], [14j . |15j . [16j . 
[17j and of H. Abidi in [1] and Abidi, Paicu in [2j. All these works give results of strong 
solutions in finite time in critical spaces for the scaling of the equations. More precisely 
R. Danchin show the existence of strong solution in finite time in [16j, [Hj, [18j when the 

initial data check (p^^ - 1, uq) G (-82% n L°°) x B.^^^ or (pg ^ - 1, uq) G Bj^-^ x B^^^^ with 
1 < p < N . Moreover R. Danchin need of a condition of smallness on the initial data, it 
means that ||po ^ ~ ^11 ^ assumed small. More recently H. Abidi and M. Paicu in [2] 

N — — 1 

have improved this result by working with initial data in B^^ ^ x B^^ ^ with pi and p2 
good choosen and checking some relations. In particular they show that we have strong 

solution for initial data uq in B^^ ^ with 1 < P2 < 2A^ which improves the results of 
R. Danchin. All these results use of crucial way the fact that the solution are Lipschitz. 
In particular, it explains the choice of the third index for the Besov space where r = 1, 

JV 

indeed it allows a control of \7u in L}j^{B^^) which is embedded in L}j.{L°°). This control 
is imperative in these works to get estimate via the transport equation on the density. 

JV JV 

However the scaling of (II. 3p suggests choosing initial data (po, uq) in B^^ ^1 xBp^^r (see the 

definition of Besov spaces in the section [2]) with (^1,^2) G [l)+oo)^ and (r,r') G [l,+oo]^. 
Indeed it seems that it is not necessary by the study of the scaling of the equation to 
impos that r, r = 1 as in the works of H. Abidi, R. Danchin and M. Paicu. The goal of 
this article is to reach the critical case with a general third index for the Besov spaces 
r and r . More precisely in the sequel we will restrict our study to the case where the 
initial data {po,uo) and external force / are such that, for some positive constant p: 

— +£ — ~ — — 1 

(po -p)£ 5p7,oo n L°°, G B;ir and / G LI^{M+ , G B^lr ). 

with r G [l,+oo], e > and pi, p2 will verify appropriate inequalities ( for a definition 
of see section [2]). 

In this article we improve the result of H. Abidi, R. Danchin and M. Paicu by working 

JV 

with initial data in -BpaV with r G [1, +00]. In particular we generalize to the case of the 
Navier-Stokes incompressible dependent density the result of Cannone-Meyer-Planchon 
in [7]. For making we introduce a new idea to control the density via the transport 
equation when the velocity is not Lipschitz. We use some new a priori estimates on the 
transport equation when the velocity is only log Lipschitz. The difficulty is to deal with 
the loss of regularity on the density, that is why we need to ask more regularity on po 
to compense this loss of regularity. The crucial point is that the density stay in a good 
multiplier space. Moreover we improve the results of H. Abidi, R. Danchin and M. Paicu 
because we do not need to assume some condition of smallness on the initial density. In 
[TT] , [12] , [13] , [H] , [IS] , [IS] 5 [IZ] , PP and [2] , we need to make the additional assumption 

JV 

that p — p \s small in B^^. The reason why is that we handled the elliptic operator in 
the momentum equation of p.ip as a constant coefficient second order operator plus a 
perturbation introduced hj p — p which, if sufficiently small, may be treated as a harmless 
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source term. For smoother data however, the additional regularity can compensate large 
perturbations of the constant coefficient operator. This fact has been used in leads 
to local well-posedness results for smooth enough data with density bounded away from 
zero. The price to pay however is that assuming extra smoothness precludes from using 
a critical framework. In fact to treat the case of large initial data on the density, we 
follow a idea of R. Danchin in [?] use for the case of Navier-Stokes compressible where to 
get good estimate on the elliptic operator in the momentum equation, we need to split 
the behavior high frequencies and low-middle frequencies of the viscosity terms. 
In the present paper, we address the question of local well-posedness in the critical 
functional framework under the sole assumption that the initial density is bounded away 
from and tends to some positive constant p at infinity (or has average p if we consider 
the case of periodic boundary conditions). To simplify the notation, we assume from now 
on that p = 1. Hence as long as p does not vanish, the equations for (a = p~^ — l,u) 
read: 

dta + u • Va = 0, 

dtu + u-Vu + {l + a){VU- pAu) = f, 
divti = 0, 

ia,u)/t=Q = (ao,Mo)- 

One can now state our main result. The first theorem generalize the work of Can- 
none, Meyer, Planchon in \7\ when the third index of the Besov space for the initial data 
is in [1, +oo[. 

Theorem 1.1 Let 1 < r < oo, 1 < pi < oo, 1 < p2 < oo and e > such that: 

N N N N 

\-e < hi and 1 < — . 

Pi P2 P2 Pi 



N 



-1 ~ ^-1 

Assume that uq G Bp^^r with divuo = 0, f £ Lf^JW^ ,Bp'^.r ) and ao G Bp^^oo H L°°, 
with 1 + ao bounded away from zero and it exists c > such that: 



|oo|| n. ^ c. 



^■^ Pi ~^ P2 ^ there exists a positive time T such that system has a solution (a, u) 
with 1 + a bounded away from zero and: 



N , e 



P2 

P2,r 



a G C( [0, T] , B:^lJ ) , n G (C( [0, T] ; ) n (0, T, B^lr )) 

and Vn G L^(0,r,B 
This solution is unique when jj < + 

Remark 1 This theorem can be applied for non degenerate viscosity, but for the sim- 
plicity in the sequel the viscosity is assumed to be a constant p > 0. 

Remark 2 As in the work of H. Abidi and M. Paicu in ^2], we are able to get strong 
solution when uq G -Bpl.r with 1 < P2 ^ 2N , it improves the result of R. Danchin in 
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USI and 173/. 

Moreover we get weak solution with initial data very close from (oo,uo) G B'n'^ ^ ^ooV 
and (ao,iio) £ -^00,00 ^ ^Nr- means that in the first case we are not very far from the 
Koch-Tataru in 126^ initial data for the velocity uq, and in the second case we are close 
to just impose a condition of norm L°° on oq which is of great interest for the system 
bifluid. 

Remark 3 In the previous theorem, we need of a condition of smallness, because when 
P2 / 2, we have extra term in our proposition [STSl which requires a condition of smallness 
on a. It seems that there is a specific structure in the scalar case when p2 = 2. 

In the following theorem, we improve the previous result in the specific case where p2 = 2. 
In this case we don't need to impose condition of smallness on the initial data. 

Theorem 1.2 Let 1 < r < 00, 1 < pi < 00 and e > such that: 

N N N N 

— +£< — + 1 and —<! + —. 
pi 2 2 pi 

K-i ~ ^-1 —+e 

Assume that uq e with divno = 0, / G Lj^^{R'^ , B^^ ) and ao G -Bpi,oo n L°°, 

with 1 + oo bounded away from zero. There exists a positive time T such that system 
has a solution (a, u) with 1 + a bounded away from zero and: 



N 



JV_i ^ M 



and Vn G L'{0,T,B. 



N__ 

l/n R 2 

2,r 



1, 



This solution is unique when < + ^. 

In the following theorem we choose r = +00, we have to treat the case of a linear loss of 
regularity on the density p which depends of the solution u. This fact come from your 

^ N_ 

estimate on the transport equation when Vu G L\,[Bp^^oo), in this case we have a loss of 
regularity on the density which depends of this quantity. 

Theorem 1.3 Let 1 <pi < oo, 1 < p2 < 00, and e > such that: 

N N N N 

he < hi and — < H . 

Pi P2 P2 Pi 

^-1 , ^-1 

Assume that uq G Bp^^oo with divwo = 0, / G L;^^(M+, i3p|^oo ) and oq G Bp^,oo H L°°, 
with 1 + ao bounded away from zero and it exists c > such that: 

||oo|| M. < C. 

,00 

// ^ + ^ > j^, there exists a positive time T such that system {T^P ^ solution (a, u) 
with 1 + a hounded away from zero and: 



G C7([0, T], u G (C([0, T\-B;Ioo n Li(0, T, B^.^oo )) 



N 



and vn G L^{{),T,B^^ 



^~i 
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with: 



o-(r) = — + e- A||n||^ jv+i +A / W{t)dt, 

L\.{B^loo ) Jo 



for any A > and any nonegative integrahle function W over [0, T] such that cr{T) > 
— 1 — A^min(^, 4-). This solution is unique when jj' < ^ + 

In the following theorem, we generalize the previous result with large initial data for the 
density when p2 = 2. 

Theorem 1.4 Let 1 < pi < cxd and e > such that: 

N N N N 

— +£< — + ! and —<! + —. 
Pi 2 2 pi 

Assume that uq G Ti2oo with divuo = 0, / G L\^J^'^ ^B^ ^ ) and ao G -Bpi,oo H L°°, 
with 1 + oo bounded away from zero. There exists a positive time T such that system 
has a solution (a, u) with 1 + a bounded away from zero and: 



1, 



, , 1\_ _ -| JV_ 1-1 J, J 

aGC([0,r],S,-(S), nG(C7([0,r];i?2% n l1(0, T, + ))^ 

and Vn G L^{{),T,B^^^ 
with: 

o-(T) = — + e - AllnlL n +X W(t)dt , 

Pi L},{B^^') Jo 

for any A > and any nonegative integrable function W over [0, T] such that (t{T) > 
— 1 — iVmin(^, \). This solution is unique when < ^ + \- 

The key of the theorems 11.11 \\.2\ 11.31 and 11.41 is new estimate for transport equation for 
u which is not considered Lipschitz. In this case we have to pay a loss of regularity on 
the density, p. The basic idea is to deal with this loss of regularity by conserving a 

in CxiBp^ _^^) n L°° which have good properties of multiplier. Moreover to avoid as in 
the work s of H. Abidi, R. Danchin and M. Paicu conditions of smallness on the initial 

density, the basic idea is that having coefficients in a in CxiB^^^j^^) provides us with 
some uniform decay on the high frequencies of the variable coefficients so that the elliptic 
operator may be considered as an operator with smooth coefficients (of the type which 
has been investigated in [15]) plus a small error term. 

Our paper is structured as follows. In the section [2l we give a few notation and briefly 
introduce the basic Fourier analysis techniques needed to prove our result. Section [3] and 
Hlare devoted to the proof of key estimates for the linearized system (II. 4p in particular the 
elliptic operator of the momentum equation with variable coefficients and the transport 
equation when the velocity is not assumed Lipschitz. In section[5l we prove the existenceof 
solution for theorem ll.ll and ll.2l whereas section [6] is devoted to the proof of uniqueness. 
Finally in section [71 we briefly show how to prove theorem 11.31 and 11.41 Elliptic and 
technical estimates commutator are postponed in an appendix. 
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2 Littlewood-Paley theory and Besov spaces 



Throughout the paper, C stands for a constant whose exact meaning depends on the 
context. The notation A < B means that A < CB. For all Banach space X, we denote 
by C([0, T],X) the set of continuous functions on [0, T] with values in X. For p G [1, +oo], 
the notation LP{0,T,X) or Lj,{X) stands for the set of measurable functions on (0, T) 
with values in X such that t \\fit)\\x belongs to LP{0,T). 



2.1 Littlewood-Paley decomposition 

Littlewood-Paley decomposition corresponds to a dyadic decomposition of the space in 
Fourier variables. Let a > 1 and {(p,x) t»c a couple of smooth functions valued in [0, 1], 
such that ifi is supported in the shell supported in G M^/a"^ < |^| < 2a}, x is 
supported in the ball G < a} such that: 

Denoting h = T'^ip, we then define the dyadic blocks by: 
Aiu = if / < -2, 

A_in = x{D)u = h* u with h = ^~^x, 

Aiu = ip{2-^D)u = 2'^ / h{2^y)u{x - y)dy with h = T'^x, if / > 0, 
Siu = ^ Afcit . 

h<l-\ 

Formally, one can write that: u = X^^g^ ^fe"" • '^'^^^ decomposition is called nonhomoge- 
neous Littlewood-Paley decomposition. 



2.2 Nonhomogeneous Besov spaces and first properties 
Definition 2.1 For s G M, p G [1, +oo], q G [1, +oo], and u G >S'(M^) we set: 

lei 

The Besov space B^^ is the set of temperate distribution u such that ||n||Bs^ < +oo. 

Remark 4 The above definition is a natural generalization of the nonhomogeneous Sobolev 
and Rdlder spaces: one can show that B^^^ is the nonhomogeneous Rdlder space and 
that B2 2 is the nonhomogeneous space . 

Proposition 2.1 The following properties holds: 

1. there exists a constant universal C such that: 
C~^\\u\\b'> < I|V«|| 753-1 < C||n||Rs . 

II II p,r — II ii-Op^T- — II ll-^p,r 

2. Ifpi < p2 and n < r2 then B;^^,^ ^ B^-p^^'^'-^'^'l 
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Let now recall a few product laws in Besov spaces coming directly from the paradifferen- 
tial calculus of J-M. Bony (see |6]) and rewrite on a generalized form in [2] by H. Abidi 
and M. Paicu (in this article the results are written in the case of homogeneous sapces 
but it can easily generalize for the nonhomogeneous Besov spaces). 

Proposition 2.2 We have the following laws of product: 

• For all s G M, (p, r) S [1,+cx)]^ we have: 

\\uv\\^s^ <C{\\u\\loo\\v\\bs^^ + \\v\\l^\\u\\b;^^^). (2.5) 

• Let {p,pi,P2,r,Xi,X2) G [1, +00]^ such that:^ < ^ + j^, Pi < X2, P2 < Ai, i < 
^ + jj- and p ^ ^ + ij- ^6 have then the following inequalities: 

^fsl + S2 + A^inf(0, 1 - ± - ±) > 0, si + < ^ and S2 + < then: 

||n7;|| < Mb;iJv\\b;i^, (2-6) 

when si + g = ^ (resp s2 + ^ = ^)we replace \H\b;iM\b;1^ T^sp \\v\\j^p^ J 

by \W\\B;iJHB;ir T^sp \HB;i^nL^)' + S = f ^2 + = we take 
r = l. 

// si + S2 = 0, si e (f - ^, ^ - f 1 and^ + ^ <1 then: 

_iv(i + i-l) ^ Ihllfi^i Iblls^a . (2.7) 

// \s\ < ^ for p > 2 and —y < s < ^ else, we have: 



\uv\\b^^,. < C\\u\\bs J\v\\ n . (2.8) 



Remark 5 In the sequel p will be either pi or p2 and in this case j = ~ ^ "if Pi ^P2, 

Corollary 1 Let r G [1, +00], 1 < p < pi < +00 and s such that: 

. s G (-^,^) + ^ < 1, 
\ Pi ' Pi ' '' p Pi — ' 

• «e(-f + ^^(i + ^-l),f)^/i + J^>l, 

JV 

then we have if u £ Bp ^. and v G Bp^^oo H L'^ : 

WuvWb" < C'll'u||_Rs \\v\\ N 

II II ^p,r — II II -^p,r II II 



8 



The study of non stationary PDE's requires space of type LP{0,T, X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to 
localize the equation through Littlewood-Payley decomposition. But, in doing so, we ob- 
tain bounds in spaces which are not type LP{0, T, X) (except if r = p). We are now going 
to define the spaces of Chemin-Lerner in which we will work, which are a refinement of the 
spaces Lj,{Bp ,^). 



Definition 2.2 Let p E [l,+oo], T G [l,+oo] and si £ R. We set: 



We then define the space Lj.(Bp^j.) as the set of temperate distribution u over (0, T) x 
such that \\u\\rpfn=i < +00. 

We set CriB^^r) = ^t(-^p,V) n C([0, T], S^i^). Let us emphasize that, according to 
Minkowski inequality, we have: 



Remark 6 It is easy to generalize proposition \2.SX to Lf^{Bp^^) spaces. The indices si, 
p, r behave just as in the stationary case whereas the time exponent p behaves according 
to Holder inequality. 

Here we recall a result of interpolation which explains the link of the space B^ ^ with the 
space Sp oo) see [TT] . 

Proposition 2.3 There exists a constant C such that for aZZ s € M, e > and 1 < p < 

+00, 

II II ^ n^ + ^w II All ^^^t{^^\ 

\\u\\lP,Bs ^<C \\u\\lP,jjs Jl+logf—f^ . 



Definition 2.3 Let F be an increasing function on [1,-|-cxd[. We denote by Br{R^) the 
set of bounded real valued functions u over MJ^ such that: 

||VS','u||ioo 

\mBr = \\u\\l°° +SUP z:^— T < +00. 

We give here a proposition concerning these spaces showed by J-Y. Chemin, see [5]. 

Proposition 2.4 Let e > and u E L\.{Bp^r ) then we have u G L\,{By{M.^)) with 
T{t) = {-\ogtY+^--r forO<t<l . 
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3 Estimates for parabolic system with variable coefficients 



In this section, the fohowing linearization of the momentum equation is studied: 

f dtu + 6(Vn - /iAu) = f + g, 



divM = 
L u/t=o 



(3.9) 



Uq. 



where b, f, g and uq are given. Above u is the unknown function. We assume that 
Uq G Bp ,^ and / G -^^^(0, T; Bp ^.), that b is bounded by below by a positive constant b and 



that a = b — 1 belongs to L°°(0, T; Sp/^oo ) H L°°. In the present subsection, we aim at 



proving a priori estimates for ()3.9p in the framework of nonhomogeneous Besov spaces. 
Before stating our results let us introduce the following notation: 



At =l + b-^\\yb\\^ jv+,_i witha>0. 



(3.10) 



Proposition 3.5 Let !^ = bp, and {p,pi) G [l,+oo]. 

. Ifpi >pwe assume that s G (-^, f ) + ^ < I and s e (-^+Ar(i + J-_i), K) 
if- + ->l. 

• -^f Pi ^ P ^/len we suppose that s G (— '^) if P > 2 anc? s G (~^; ^/P < 2. 
If p ^ 2 we need to assume than there exists c > such that: 

||a|L^ <c. 

Lei m ^ Z be such that 6m = 1 + •S'mO satisfies: 

inf 6m(t,2;)>|. (3.11) 

There exist three constants c, C and k (with c, C, depending only on N and on s, and k 
universal) such that if in addition we have: 

||« — 'S'mfllL ^ ^ c= (3.12) 

then setting: 

Zrait) = 2^"'"fl\-^ [ \\af jv dT, 

Jo B^l^rM°- 

Let a' > checking d < min(l,a, ) . We have for all t G [0,T] and k = A-: 

ML^iB^p^) + '^J^ML},iB;t')^^''^^^^^ 

E.{p-1). 



i_ 

TL 

Z-(b;:/"-) - ' p ^" "L\{Biy-'^ y 



+ ^^^\\Vg\\^._,+^+^^^{ — - — )^t||u|| 
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Moreover we have VII = VIIi + with: 



6||vn2||n 



<AU\\Qg\L 



T m ( ft 



2_ + u a 



(-Bp 



^ )• 



(3.14) 



^ N_ 

Remark 7 Xei us stress the fact that if a £ L°°((0, T) x Bpl^oo) then assumption ^3.11\) 
and i3.12\) are satisfied for m large enough. This will be used in the proof of theorem 
and \1.4\ Indeed, according to Bernstein inequality for m large enough S3.11\) and ^3.12^) 
are satisfied. 

Proving proposition 13.51 in the case b = cste is not too involved as one can easily get rid 
of the pressure by taking advantage of the Leray projector V on solenoidal vector-fields. 
Then system ()3.9p reduce to a linear ipDO which may be easily solved by mean of energy 
estimates. In our case where b is not assumed to be a constant, getting rid of the pressure 
will still be an appropriate strategy. This may be achieved by applying the operator div 
to (|3.9p . Indeed by doing so, we see that the pressure solves the elliptic equation: 

div(6Vn) = divF. (3.15) 

with F = f + g + fiaAu. Therefore denoting by TCb the linear operator F — > VII, system 
(13. 9p reduces to a linear ODE in Banach spaces. Actually, due to the consideration of 
two forcing terms / and g with different regularities, the pressure has to be split into two 
parts, namely 11 = Hi + 112 with: 

div(6ni) = div/ (3.16) 

div (6112) = divH and H = g + fiaAu. (3.17) 
Proof of proposition 13. 5t 



Let us first rewrite ()3.9p as follows: 

' dtu - bmfiAu + bVU = f + g + Em, 

< divn = 0, (3.18) 

^ Ut=0 = Uq. 

with Em = fj.Au (Id — Sm)ci and 6m = 1 + SmO-- Note that by using corollary [1] and as 
— ^<s<;^forp>2or4-<s<:^ else, the error term Em may be estimated by: 

ll^mlls^^,, < ||a - 5'ma|| 2L II-D^uIIbs^^. (3.19) 

Now applying operator and next operator of free divergence yield V to momentum 
equation (|3.18p yields: 

- ;udiv(6„Vn,) = Vfg + Vgg + AgVEm + Rg - AgV{aVU), (3.20) 
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where we denote by Uq = AgU and with: 

Rq = R\+ Rq- 

where: _ 

R] = is{VAg{b^Au) - Vdw{bmVug)), 

Rq = A'('^div(6mVtig) — div(6mVuq)) = —fj,Qdiv{SmaVug). 

where Q is the gradient yield projector. 

case p ^ 2 

Next multiplying both sides by \ug\P~'^Ug, and integrating by parts in the second, third 
and last term in the left-hand side, we get by using Bony decomposition (for the notation 
see HU): 

~^ll"gllLP + / hm\'^Ug\'^\Ugf~'^dx + [ hm\ug\^~^\V\u\^\'^dx 

P dt JjRJV J^N 

< hqfLp'iWrfghp + WVgqhv + \\Rq\\LP + 1 1 A, (TvaH) 1 1 + 2^Ag{TaU)\\Lr> ^^''^^^ 

+ ||Ag(rvna)||LP + \\VAgErr,\\Lp)- 

Indeed we have as divu = and by using Bony's decomposition and by performing an 
integration by parts: 

/ Ag{aVU)\Ugf-\g = [ Ag{T^^a)\Ug\P-\dX - [ Ag{T^a'n)\Ug\P-^UgdX 

JRN Jrn Jrjv 

- / Ag{TaU)dW{\Ug\P-\)dX. 

Next we have: 

V{\Ugr^) ■ Ug = {p-2)\Ugr^Y.<^^'4<^ 

i,k 

and by Holder's and Berstein's inequalities: 

l|V(|u,r-2) . ug\\^_^ < dp - 2)2^Ugf-\ 

Next from inequality (j3.2ip . we get by using lemma A5 in [16j : 

^^ll^^lli. + ^^^^2^^Ug\\l, < \\Ugrj-'{\\Vfg\\L. + ^ <? J | + \\VAgE^\\L. 

+ ||A,(rvan)||iP + 2''||A,(r,n)||iP + ||A,(rvna)||LP + WRqhp), 

Therefore, elementary computation yield (at least formally): 

" Jt^""^ <\\'Pf,\\LP + \\'P9,\\lv + \\VAgEm\\Lv 



||A,(Tvan)|Up + 2''||Ag(T„n)||iP + \\Ag{Ty^a)\\Lv + \\Rg\\Lv. 
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We thus have: 

hMlL^ < e-^'''*||A,no||Lp + f e-'^^"'"^'-^\\\VU\\Lv + 11^5,11^.+ 

JO 

\\V\Em\\Lv + \W{TyJi)\\Lv + 2^||Ag(r,n)||iP + ||Ag(Tvna)||LP + ||ii,||Lp)(r)dT, 

which leads for ah q > —1, after performing a time integration and using convolution 
inequalities to: 

^ ~^^p2 ^^ )^'^'^\\^q\\Lr^{LP) < ||A,Uo||lp + ||7^/g||Li,{LP) + l|Ag(rvani)||^i^(^p) 

+ 2'?||A,(r„ni)||i^(i,) + ||A,(rvn,a)|L^(ip) + + ||PAgi?„||^^(^,) 



+ ll^g(^vn2«)llLi,(LP) + WPaqhipiLP)) 



We are now interested by treating the commutator term , we have then by using lemma 
[T]in the appendix the following estimates with a < 1: 

mlhr < c,P2(-i+")«^||S^a|| ^^JDuWb^^^^, (3.23) 

,oo 

where (cg)ggz is a positive sequence such that Cq £ F , and i? = fi. Note that, owing to 
Bernstein inequality, we have: 



||5'ma|| 



N 

— +a r-^ ,, ,, 



Next we have by corollary [T} 

WRIWlp <Cgi^2-'^'\\Sma\\ jv ||n||^.+2, (3.24) 



Hence, plugging (|3.23p . (|3.24p and (|3.19p in (j3.22p . then multiplying by 2'^* and summing 
up on (7 G Z in F, we discover that, for all t £ [0, T]: 

\W\\l°°(B'> ) + (— — -^)'^II^IU , < ||^^o||b» + ll^/llrwR. -) 

+ ||T„ni||2;i^(^s+i) + ll^vanillji^^^,^^^ + l|rvnia|lzi,(ij._j 

+ Ci^\\a- S„,a\\^ ^ 11^1111(5^+2) +2""" / ^ WII^^IIb^+i W^^^ (3.25) 

+ (^^^^)^(ll^fflU s-2+^ +l|ran2|U .-1+^ +||Tvan2|U s^2+^ 
^ p2 ^ll^„^^s^^2+„^ 

for a constant C depending only on N and s. With our assumption on a, a and s, the 
terms II^VaHi and l|rvni"llLi,(B|J bounded by: 

l|a|| Jv+„ llVnilL, 
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whereas llTvallolL , 9, 2 and llTvyn « 94.2 may be bounded by: 

Moreover we control llTalli llri ,„s+l^ and ||ran2|L ,_i + a by respectively: 



livniii^i 



||o|| jv ||Vn2|L ,_24--2- • 

Hence in view of proposition 18.81 and provided that < a' < min(l,a, |) and < ^ 
(which is assumed in the statement of proposition 13.50 and a G [0, a ], 



^ll^n^llz^(Br/')--^^" (3-26) 

On the other hand, by virtue of proposition 12.21 and of assumption on a, a , a" and s, 
we have: 

L-(B;,'+- ) L^iB;,r^-- ) L^{Bp\\tnL^) LrfiBl,r^^ ) 

As a < min(l,a, |(s — 2 + ^)), proposition 18.81 with (T = s — 2 + ^— a" (here comes 
s > 2 — ^) applies, from which we get for all e > (e = does if m > 2), 

L™(Bp,, ^"^ ) L™(iJj,,, +™ (3.27) 

Let X(t) = ||tt||/,oo(gs ^-j + z^6||u||j^ij^^a+2-). Assuming that m has been chosen so large as 
to satisfy: 

CuWa — Sm^W jv < z^, 

i?f(Bjf/,oonL°°) 

and by interpolation we get: 

^2p222ma 

Cz^2™"||a|| JV ||u||r.+2 < + ||a||^ iv \\u\\b^ , (3.28) 

Plugging ([OH]) . (fOT]) and (lOHl) in ([3:25]) . we end up with: 

X{T) < Ikobi,. (r/|lzj(^.„) + + (f 2^ 

X X(T))dr + ( -^^ ~ )^^r||n||-, ^2 'Y 
Gronwall lemma then leads to the desired inequality. 



^^™"'|af JV (r) 

-'^P]^ ,00 
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Case p = 2 

In this case we don't need of condition of smallness on ||a|| jv , indeed the bad 

R PI n,T oc 
^Pl ,ool 1-^ 

terms as or 2'?||A5(Tan)||/^2 disapear in the integration by parts as divug = 0. So we 
can follow the same procedure and conclude. □ 

4 The mass conservation equation 

4.1 Losing estimates for transport equation 

We now focus on the mass equation associated to p.ip : 

{dt.a + V ■ Va = g, ^ ^ 

(4.29) 
a/t=o = ao- 

We will precise in the sequel the regularity of oq, v and g. In this section we intend to 
recall some result on transport equation associated to vector fields which are not Lipschitz 
with respect to the space variable. Since we still have in mind to get regularity theorems, 
those vector field cannot be to rough. In order to measure precisely the regularity of the 
vector field v, we shall introduce the following notation: 

7 — 

t<..W=sup ^''l«ll^" <+oo. (4.30) 

j>0 U "T ^) 

Let us remark that if pi = +oo then Vp_^ ^ is exactly the norm ||i?r|| of definition 12. 3i 
4.1.1 Limited loss of regularity 

In this section, we make the assumption that there exists some a s]0, 1[ such that the 
function „ defined in (|4.3U|) be locally integrable. We will show that in the case a = 1, 
then a linear loss of regularity may occur. In the theorem below, Bahouri, Chemin and 
Danchin show in [5j that if a g]0, 1[ then the loss of regularity in the estimate is arbitrarily 
small. 

Theorem 4.5 Let {p,pi) be in [l,+oo]^ such that 1 < p < pi and a satisfying a > 
-1 - iVmin(^,^). Assume that a < 1 + ^ and that T^^ ^ e]0, l[ is in L^{[0,T]). 

Let ao G Bp,^ and g £ L\,{Bp^^). Then the equation ^.29\ ) has a unique solution 



a E C([0, T], n^/^^i?p and the following estimate holds for all small enough e: 

W^h^iB^r^) - ^(ll«0b-^ + hWl^^^B'^^^)) exp (^(^pi,a(T))^), 

where C depends only on a, p, pi, a and N . 

In the following proposition, we are interested in showing a control of the high frequencies 
on the density when u is not Lipschitz. Indeed we recall that in the proposition 13.51 when 
p = 2, we need to control the high frequencies of the density. In particular the following 
proposition is useful only in the case of theorem 11.21 
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Proposition 4.6 Let {p,pi) be in [l,+oo]^ such that 1 < p < pi and a satisfying 
a > -1 - iVminf^, \). Assume that a < 1 + — and that V' „ gIO, 1[ is in L^ffO-Tl). 

Let aQ £ and g G L}j.{Bp^^), the equation \4.29 ) has a unique solution a £ 



C{[0,T],n^i^^Bp^) and the following estimate holds for all small enough e: 
5]2('^-)'||A,a(t')|U^(iP)<^(2'^'||A,ao||Lp) + Cr?i^^V,;,„(t') 



l>m l>m 



X (Ikolls-^ + ||s'|l2;i,(sa^)) exp {^^{Vp^,a{t))^)dt . 



C 

e 

where C depends only on a, p, pi, a and N . 
Proof: 

By using the proof of Bahouri, Chemin and Danchin in 5j one can write: 

2(2+^>*||A,a(0||L. <2(2+0<x||A^;^||^^ + C7(— — / Vp^^J)\\a{t )\\ .,dt. (4.31) 
Whence taking the supremum over I > m, we get 

sup sup(2'^*''||A,a(t')||Lp) <sup(2'^'||Aiao||Lp) + C7r/i^ / Vp^^^{t)\\a{t)\\ ydt. 

We apply Gronwah inequality in ()4.3ip and we insert in previous inequality: 

sup sup(2V'||A;a(t')||Lp) < sup(2'^' || A^aoHip) + Ct?^ ||/o|bj^,^ t K.A^') 

t' G[0,t] l>m Jo 

C / 1 / 

X exp (^^iVp-,.ait))^)dt . 

which leads to the proposition. □ 
Remark 8 In the sequel, we will use the theorem |^.5| and the proposition \4.6\ when 

pi = oo and a = 1 + e — - ■ Indeed we will have u is in L {Bp^^r ) (with 1 < r < cc) 
and according proposition \2.4\ and definition \2.3\ we get: 

ft 



VooAt )dt < ||n|L 



JV 



-1 



So it will allow to get estimates on the density with an arbitrarly small loss of regularity. 
4.1.2 Linear loss of regularity in Besov spaces 

Tis section is devoted to the estimates with linear loss of regularity. Remind that u is 
log-lipschitz only if only if there some constant C such that: 

||V5j||loc < C{j + 1) for all j > -1. 

So Bahouri, Chemin, Danchin show in ^ the following theorem. 
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Theorem 4.6 Let 1 < p < pi < co and si G M satisfies a > — A^min(^, \). Let a in 

]si, 1 + ^[ and V a vector field. There exists a constant C depending only on p, pi, a, si 
and N such that for any X > C, T > and any nonnegative integrable function W over 
[0,T] such that if ax > si with: 

at = a-X f [V^,,i{t) + W{t))dt 
Jo 

then the following property holds true. 

Let oo S -Bp oo and g = gi + g2 with, for all t G [0, T], gi{t) G -Bp'oo ^''^d: 



Vj > -1, ||A,52(t)||LP < 2-^'^^{2+j)Wit)\\fit)\\^^^.^. 
Let a G C{[0,T]; Bp^^) be a solution of \4.29\) . Then the following estimate holds: 

A '■^ 

ie[o,T] 



sup llaWlb-t^ < Y^{\\ao\\B-^^ + 



I \\gi{t)\\B^^^^dt). 
Jo 



In the following proposition, we generalize this result to the high frequencies of the density 
P- 

Proposition 4.7 Let 1 < p < pi < oo and si G M satisfies a > — A^min(^, jr). Let a 



in ]si,l + — [ and v a vector field. There exists a constant C depending only on p, pi, 



-I 

Pi I 

a, si and N such that for any X > C, T > and any nonnegative integrable function W 
over [0, T] such that if gt > si with: 

ft 



at = a-X f {Vp^^{t) + W{t))dt 
Jo 



then the following property holds true. 

Let oq G -Bp^oo ^''^d g = gi + g2 with, for all t G [0, T], gi{t) G i?p*oo 

Vj > -1, ||A,52(t)||LP < 2-^'^\2+ 3)W{t)\\fmB;u- 
Let a G C([0, T]; 5^^^^) be a solution of l^4-^^ )- Then the following estimate holds: 



X r'^ 

sup ^2''^iA/a||LP < — ^(llaob.^ + / \\giit)\\B^ dt). 
t&,Tf^ X-C 'Jo 



Proof: 



The proof follows the same line as in proposition 14.61 □ 
Remark 9 In the sequel, we will use the theorem \4-&\ cmd the proposition \4. 7| when 



pi = oo. Indeed we will have u in L [Bp^^oo ), then according proposition \2.1l Vu is in 
L^{B^ ^)so we obtain the following control: 

SUp( p— )dt < IL+i ■ 

l>0 i + 1 L^B^l^) 

So it will allow to get estimates on the density with a small loss of regularity on a small 
time interval. 
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5 Proof of the existence for theorem 11.11 

We use a standard scheme: 

1. We smooth out the data and get a sequence of smooth solutions {a'^,u^)n£N to 
(II. 4p on a bounded interval [0, T"] which may depend on n. 

2. We exhibit a positive lower bound T for T", and prove uniform estimates in the 
space: 

Et = Ct{B^IJ) X {CriB^lr ) n L^(S/„. )))''. 
for the smooth solution {a^^u^). 

3. We use compactness to prove that the sequence converges, up to extraction, to a 
solution of (fT7 



Construction of approximate solutions 

We smooth on the data as follows: 

Oq = S'nfflo, = SnUQ and /" = Snf. 

Note that we have: 

V/ e Z, IIA^OoIIlpi < IIA/OoIIlpi and ||ao|| n_^^i < \\ao 



and similar properties for Uq and a fact which will be used repeatedly during the 
next steps. Now, according [2], one can solve p.4p with the smooth data {oqjUq,/^). 
We get a solution (a",ii") on a non trivial time interval [0,T„] such that: 

(b.62) 

and Vn" G L^^iH^). 

Uniform bounds 

Let Tn be the lifespan of (a„, u„), that is the supremum of all T > such that (jl.4p with 
initial data {u^jUq) has a solution which satisfies ()5.32p . Let T be in (0,T„), we aim at 
getting uniform estimates in Et for T small enough. For that, we need to introduce the 
solution (uj, Vn^) to the nonstationary Stokes system: 



'dtul-fxAul + VUl = r, 

dWul = 0, 
. K)t=o = 



Now, the vectorfields vP' = vP" — and VH" = VH^ + VII" satisfy the parabolic system: 

' dtu"" - ^(1 + a")An" + (1 + a")Vn" = F(a", n", VH"), 
< divn" = 0, (5.33) 
. ^i"(0) = 0, 
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which has been studied in proposition 13.51 with: 

//(a", u", Vn") = a"(/uA< - VHJ) - • Vw". 

Define m G Z by: 

m = mf{peZ/ 2iy'^2^n\\Aiao\\LPi <cu} (5.34) 



l>q 

where c is small enough positive constant (depending only N) to be fixed hereafter. 
Let: 



6 = 1+ sup ao(x), Aq = 1 + 2\\ao\\ n_ ' , Uq = \\uo\\ n_.i + 

and Uq = 2CUq + ACvAq (where C stands for a large enough constant depending only 
N which will be determined when applying proposition 12.21 13.51 and 14.61 in the following 
computations.) We assume that the following inequalities are fulfilled for some r/ > 0, 
a > 0: 

{Hi) ||a"-5ma"|L jv <cvp-^, 



(Wa) Cp2r||a"f JV <2-2"^z/, 

-'T \ Pi ' 



(W4) ||a"L^ ^ < Ao, 

(W5) IKIL i^+i 

(^e) ||S"|L iv„i +id||5"IL < f/or?, 

(W7) l|vn2IL iv_, <r?, 

(Ws) l|vn"|| <no??, 

^P2. 



Remark that since: 

1 + Sma"" = 1 + a" + (5^a" - a"), 
assumptions (Wi) and (Ws) insure that: 

inf (l + 5™a")(t,x) > ^6, (5.35) 
provided c has been chosen small enough (note that | < 6). 

We are going to prove that under suitable assumptions on T and rj (to be specified 
below) if condition ijii) to iJ~L%) are satisfied, then they are actually satisfied with strict 
inequalities. Since all those conditions depend continuously on the time variable and 
are strictly satisfied initially, a basic boobstrap argument insures that iTLi) to ("Hs) are 
indeed satisfied for a small r„ < T". First we shall assume that r] satisfies: 

C {I + uT^Uq)!] <\og2 (5.36) 
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so that denoting: 

ft 



(V^")p„i_i(t)= iiVn^^^^_,{s)ds and = L iK\,,,.i{s)ds, 

with: 

l>0 + 

We recall now that according proposition 12.41 

iVnr,„i-dt)<C\\u^__ and (FfL . i (t) < C||n2|L 



we have, according to (Hs) and (We): 



(5.37) 



In order to bound a" in L'^{Bpl^oo )■, we apply theorem 14.51 and get: 

ll«"ll Jv,. < e^(i)'"7||n"|| jv +IKII jv )||ao|| jv,,. (.38) 

Moreover as we know that ||a"||ioo < ||ao||L°°, is satisfied with a strict inequality. 
Next, applying classical proposition on heat equation ( see [H]) yields: 

"V jv_i <t/o, (5.39) 

i^rp \ijp2 ,r } 

N 



^Hul\\_^ < (^2'^^^i-^^(l-e— 2 ^)'-(||A,no|ri.2 +||Az/||2i(M+LP2)))". 

(5.40) 

Hence taking T such that: 

( j;2''^(^-^)(l - e~^^^''^ymM\h + l|A//|l2i(M+,L^)))' < ^r,y, (5-41) 



insures that (7^5 ) is strictly verified. 

Since (Wi), ("^2) and ()5.35p are satisfied, proposition 13.51 may be applied with a 
We get: 



IS" I 



|a"(A<- vn2)|L +||7x"-Vu"|L jv_i +^T,n||'u|L 



20 



with Z^{T) = T^^v^v'^ l|a"f jv cLt. Next using Bony's decomposition and 
divn" = 0, one can write: 

with the summation convention over repeated indices. 

Hence combining proposition 1.4.1 and 1.4.2 in [llj with the fact that L^^Bp^^r ^) ^ 
~ _i 

L^(Soofoo) for p = I or p = 4, we get: 

||div(u"8)n")|L <C||n'^|L4 i^+i lh"IL f-i- 

By taking advantage of proposition 13.51 12.21 12.11 and Young'inequahty, we end up with: 



||«^IL +k\\u''L^ < e^^^^^^ (\\uo\\ i|_i + (||a"|L + 1)^ 

X (CIKIL ^+1 (z^lla^ll JV +IKIL jv_j + ||vn2iL jv.^ x 

||a"|| |t ) + (l|a"|L l^+f +1)^*I|2"IL^ n+i]. 

with C = C{N) and Cg = {N,g,b,b). Now, using assumptions (7^2 )i ("^4), ("^5)1 ("^e) 
and (Ti-j), and inserting ()5.37p in the previous inequahty and choosing T enough small 
gives: 

||S"|L +I|S'^IL f+i <2C{i?Ao + Uo)r] + 2CyTAo, 

Lrp (^Bp2^r ) Lrj-i(^Bp2,r ) 

hence (We) is satisfied with a strict inequality provided: 

C7gr < CvT]. (5.42) 

To show that {Tij) and ("Hs) are strictly verified on (0,T^), we proceed similarly as for 
(7^5) and iJ~L&). We now have to check whether iJ~Li) is satisfied with strict inequality. 
For that we apply proposition 14.61 which yields for all m £ 7^, 

1 N 1 N 

sup2^i||A;a"||is?(LPi) <C sup2^i||A/ao||LPi) (||^i"IL ^ + IKIL Jv (5.43) 

l>m ' l>m LliB^lr) L\{B^lr) 

Using (|5.36p and (Ws), (We), we thus get: 

lIL C _i~ 
Wa^ - SmO^W JV < SUp2 Pi ||A;ao||LPi + :j ^(l + ||«o|| JV + U '^Lq)!]. 

i§?{B^fl^oc)nL°° l>m log 2 B^l^nL^ 

Hence (TCi) is strictly satisfied provided that rj further satisfies: 

^ -Al + WaoW N )(l + z.-iLo)ry< ^. (5.44) 



log 2 B^l.^nL^ 2i/ 



So TCi is strictly verified. 

(^3) is trivially verified by the transport equation as we assume that 1 + oq is bounded 
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away and that oq G L°° ■ 

Next, according to (7^4) condition ("^2) is satisfied provide: 

One can now conclude that if T < has been choosen so that conditions (I5.4ip . ()5.42p 
and (j5.45p are satisfied (with r] verifying (|5.36p and (j5.44p . and m defined in (j5.34p and 
n > m then (a", u"", H") satisfies (TCi) to [Tig), thus is bounded independently of n on 

[o,r]. 

We still have to state that may be bounded by below by the supremum T of all times 
T such that ()5.4ip . ()5.42p and ()5.45p are satisfied. This is actually a consequence of the 
uniform bounds we have just obtained, and of a theorem of blow-up of R. Danchin in 
|14j . Indeed, by combining all these informations, one can prove that if T" < T then 
(a", u"-, Vn") is actually in: 

hence may be continued beyond T (see the remark on the lifespan in [T3] where a control 
of Vu in L-'^(i?^ is required). We thus have T'^ > T. 



Compactness arguments 

We now have to prove that (a",u"')nGN tends (up to a subsequence) to some function 
{a,u) which belongs to Et and satisfies (II. 4p . The proof is based on Ascoli's theorem 
and compact embedding for Besov spaces. As similar arguments have been employed in 
|12j or [13], we only give the outlines of the proof. 

• Convergence of {aP')n&^'- 

We use the fact that a" = a" — Oq satisfies: 

= -n" • Va". 

~ ^+1 ~ 

Since (u"')„gN is uniformly bounded in L\.(Bplr ) H L^{Bp^r ), it is, by interpo- 

lation, also bounded in Lj,{Bp2^r ) for any r G [1, -|-cxd]. By taking r = 2 and 
using the standard product laws in Besov spaces, we thus easily gather that (9to") 

is uniformly bounded in L^{Bpl^oo )■ 

\\dtS!^\\ N , < jv ||Va"|| N_, . 

t2(t>V1 T2lp>P2 \ ToalnVl n 

Hence (a"')„gN is bounded in L'^{Bpl^oo H Bpl^oo) and equicontinuous on [0,T] 

JY._i ——1 — ——1 

with values in Bpl^oo ■ Since the embedding Bpl^oo H Bpl^oo ^ Bpl^oo is (locally) 

JV 

compact, and (ao)„gN tends to ao in Bp^^oo, we conclude that (a")„gN tends (up to 
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extraction) to some distribution a. Given that (q,"^)^^ eN is bounded in L^{Bpl^r ), 

— +1 

we actually have a £ L^{Bp^^r )■ 

• Convergence of (u2)neN: 

Prom the definition of and classical proposition on Stokes equation, it is clear 
that (n2)n6N and (Vn2)neN tend to solution and Vn^, to: 

mL^{B;ir )nL^TiBplr ) for K)„gN and ^^(B;!,, ) for (VH^Wn- 

• Convergence of {u'^)n£N- 

We use the fact that: 

dtu"" = -ul ■ VS" - u"" ■ V-u" + (1 + a")A5" + a"An2 - < • V< - vfi". 
As (a"')neN is uniformly bounded in L^{Bpl^oo) and (u")neN is uniformly bounded 

~ — — 1 ~4 iV I 1 

in L^{Bp2r ) n L^{Bp^r ), it is easy to see that the last term of the right-hand 
side is uniformly bounded in L'^^Bp^^r ) and that the other terms are uniformly 

~4 iV 3 

bounded in Ls^Bp^r 

Hence (n"')ngN is bounded in L^{Bp^^r ) and equicontinuous on [0, T] with values in 

Bp2r +Bp2T ^ . This enables to conclude that (S"')neN converges (up to extraction) 

to some function ti G L!p(i?p2V ) ^ -^r(-^P2V )• 

We proceed similarly for (n2)neN and (n")„gN. By interpolating with the bounds pro- 
vided by the previous step, one obtains better results of convergence so that one can 
pass to the limit in the mass equation and in ()5.33p . Finally by setting u = u + ul and 
n = n -|- n^,, we conclude that (a, u, U) satisfies (|1.4p . 

In order to prove continuity in time for a it suffices to make use of proposition 14.51 In- 
deed, oq is in i?p/,oo n L°°, and having a G L'^{Bpl^^ ) n L°° and u G L\{Bp2^r ) insure 

^ JV 

that dta + u ■ Va belongs to L\,{Bpl^oo)- Similarly, continuity for u may be proved by 
using that uq G Bp^^r and that [dtu — i^Au) G L}p{Bp2^r )■ 



6 The proof of the uniqueness 

6.1 Uniqueness when I <P2 < 2A^, < ^ + ^ N > 3. 

In this section, we focus on the case > 3 and postpone the analysis of the other cases 
(which turns out to be critical) to the next section. Throughout the proof, we assume 
that we are given two solutions {a^,u^) and (a^,u^) of (11.40 which belongs to: 

^ JV I £ ^ — — 1 ~ — +1 

{C{[0,T];B^lJ)nL^) X {C{[0,T];B;ir ) H LHo,T; B^lr )f. 
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Let 6a = — , 5u = — v} and 5H = II'^ — 11^. The system for {5a, 5u) reads: 

{dt5a + ■ V5a = —5u ■ Vo"*^, 
StJ-u- (l + a2)(^A(5n- VOT) =F(a^n^ff). 



with: 

F(a^ n\ ff) = ■ V5u + 5u ■ Vu^ + 5a{fiAu^ - VH^). 

The function Sa may be estimated by taking advantage of proposition 14.51 with s 
^ 1 + We get for ante [0,r], 



\\6a{t)\\ < C\\6u ■ Vg^IL exp (^(y , (t))''), 



pi 



We have then by proposition 12.21 and [27 

||Mi)ll iv_i<C||HL iv ||ai^ exp(-^(||n2|L jv^JO' (6-47) 

For bounding we aim at applying proposition 13.51 to the second equation of (16.460 . 
So let us fix an integer m such that: 

1+ inf SmO^ > - and lla^ — n_ < c=. (6.48) 

Now applying proposition 13.51 with s = — 2 insures that for all time t G [0, T], we have: 

Ikll jv_2 +ki/||m|| jv + ||V(5n|| jv_2 <e'^^'"('^)x 

(||PF(a\n',tf)|L^ ^_2 +( -^^'~^^ .4T||n|L^ )• 



with = 2™/u2i/-i /g ||a(r)||2 ^ dr. 



Hence, applying proposition 12.21 corollary [T] and the fact that divdu = 0, we get as 
exemple: 



1 1 



\\6u-vu^\\ N_, <\\u'r N ,\w'r N^Aw^'^w ^-2+\\h\ ^ )• 

By the fact that |r<^ + ^, iV>3and^<i + ^ imply that: 
\\Sa{^iAui-V5U)\\^ ^_2 <||5a|| ^.^ (||A^xi|| jv_, 

+ ||vnH,. 

Lrp{_Bp2^r ) 

Now let choose Ti enough small to controll in (j6.49|) iv_^' and in (|6.5U|) 

ll^^^ill +||Vn-'-|| jv_;^ and by the fact that II ai II n_+^ < c with c small, 

rl|'RP2\ Til's ''2-1 To°('rP1^\ 

we obtain finally: 



71 (-rP2 \ roo/RPl -I Tl(nP2 \ 

(6.50) 

P2 



5a,5u,V6U)\\FT^ < cC\\{6a,5u,V6U)\\FT^, 
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with: 

^ JV ^ JV ^ N_—2 ~ ——2 

Ft = Ct{[0,T],B^Ioc) X (4(S;,V) nL??(i?;,V )) X ^(^pl.r )• 
We obtain so {6a, 6u, VdU) = on [0, Ti] for Ti enough small. By connectivity we obtain 
that ((5a, (5ti, (5Vn) = on [0,r]. This conclude this case. 

6.2 Uniqueness when:-|: = ^ + ^ or p2 = 2N or N = 2. 

The above proof fails in dimension two. One of the reasons why is that the product of 

JV_i IL_i jv_2 

functions does not map Bpl^oo x iJpaV -Sp2V but only in the larger space B2 This 

^ J^—i ^ J^— 2 ~ — 

induces us to bound ba in L^{Bpl^oo ) and 6u in L^(i?p|^oo ) H L}j.{Bp^^oo)- In fact, it 

is enough to study only the case ;|r = ^ + Indeed the other cases deduct from this 
case. If p2 = 2N then pi = ^as^<;^ + ^ and jj < + So it is a particular 
case of = ^ + For N = 2, we begin with P2 = 4 and pi = | and by embedding we 
get the result for l<pi<|,l<P2^4 and for 1 < pi < 4, 1 < p2 ^ |- 
Moreover in your case, it exists two possibilities, one when 1 < p2 < 2 and when p2 > 2. 
The first case is resolved by embedding so we have just to treat the case p2 > 2. We 
want show that {6a, 6u, VJII) E Gt where: 



Gt = Ct{[o, t],b;i^ ' ) X {l\{b;i^) n l^{b;i^ )) x l}t{b;i,. 



-^-1+5. 



In fact we proceed exactly as in the previous proof but we benefit that a is in Gt{\^, T], Bp^^oo ) 
which give sense to the product 5Au. It conclude the proof of uniqueness. 

6.3 Proof of theorem 11.21 

The proof follow strictly the same lines than the proof of theorem 11.11 except that we 
profits of the fact that in the proposition 13.51 in the case of p2 = 2 we do not need of 
conditions of smallness on the initial density ao • 

7 Proof of theorem 11.31 and 11.41 



In this case we have a control on u only in L^(i?p|^oo ), that is why to control the density 
in this case we need to use the proposition 14. 6[ The rest of the proofs is similar to proof 

m 



8 Appendix 

8.1 Elliptic estimates 

This section is devoted to the study of the elliptic equation: 

div(6Vn) = divF. (8.51) 

with b = 1 + a. 

Let us first study the stationary case where F and b are independent of the time: 
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Proposition 8.8 Let < a < 1, (p, r) G [l,+oo]^ and o" G M satisfy a < a < a + . 
Then the operator TCf, : F VII is a linear bounded operator in B^^ and the following 
estimate holds true: 

^llVnllfi.,, <^^^TO||QF||b^^,,, (8.52) 

with if a ^ 1: 

A = l + b~^\\Vb\\ N_ 



B. 



PI 



+ a-l • 



Proof: 

Let us first rewrite ()8.5ip as follows: 

div(6^Vn) = divF - Em, (8.53) 

witli Em = div {{Id — Sm)0'VlV). 
Apply Aq to (f8T53]l we get: 

div(6^VAgn) = divF, - Ag{Em) + Rg, (8.54) 

with Rg = div{bmVUg) - Aqdw{bmVU). Multiplying ([83^1) by AgU\AgU\P-^ and inte- 
grate, we gather: 



/ bm\VUg\^\Ug\P-^dX + / bm\V\Ug\^\PdX < ( | | divF, | | + \\Rg\\Lp)mg\\l-^ 

+ [ \{Id- Sm)a\\VUg\'^\Ug\P-^dx+ [ \{Id- Sm)a\N\^q\'^\^dx, 
Jr^ Jrn 



.55) 



Assuming that m has been choose so large as to satisfy: 

II Q II 

l|a - ^ma\\ ^ < -, 

an by using lemma A5 in [16| : 

S'^llnqlliP <2«||F,||LP + ||i?q||LP. 
By multipliyng by 2"^*^*^^^ and by integrating on F we get: 

^l|Vn||B.„ < ||Q/||B|,, + ||i?g||B|_,. 

The commutator may be bounded thanks to lemma [T] with 0<a<l, a = s — 1. We 
have get: 

^livniis.^^ < \\qf\\bs^ + ||V5^a|| jv^Jivn 

Therefore complex interpolation entails: 



s — a . 



Pi 



|vn|L.-. < ||vn|L=, l|vn||^o • 



Note that, owing to Bersntein inequality, we have: 



|V5„a|| iv^, <2-+°||a|| n . 



We have then: 



^livniiffi, < ||Q/||B|, + 2™+"||a|| ^ iivnii^^ iivnii^o • 

And we conclude by Young's inequality with pi = and P2 = ^ • And we recall 
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8.2 Commutator estimates 



This section is devoted to the proof of commutator estimates which have been used in 
section 2 and 3. They are based on paradifferentiel calculus, a tool introduced by J.- 
M. Bony in [6]. The basic idea of paradifferential calculus is that any product of two 
distributions u and v can be formally decomposed into: 

uv = TuV + TyU + R{u, v) = TuV + T^u 

where the paraproduct operator is defined by TuV = J2q Sq-iuAgV, the remainder oper- 
ator R, by R{u, v) = Agu{Aq-iv + AgV + Ag^iv) and T^u = T^u + R{u, v). 
Inequality ()3.23l) is a consequence of the following lemma: 



Lemma 1 Let pi £ [1, +oo], p G [1, +oo], a £ (1 — ^ k £ {1, ■ ■ ■ , N} and R, 
Ag{adkw) — dk{aAqw). There exists c = c(a, N, a) such that: 



q 



2'^''\\Rg\\LP <Ccg\\a\\ jv ,„||t(;|L<.+i-c (8.56) 



-"PI ,r 

whenever ——< a < — + a and where Cn £ U' ■ 

Pi — Pi 'i 

In the limit case a = we have for some constant C = C(a, N): 



2 '^^\\Rg\\Lp <C\\a\\ ^+il\\w\\ _jv+,_. (8.57) 

R PI R PI 

1 J^p,oo 

Proof 

The proof is based on Bony's decomposition which enables us to split Rg into: 

Rq = dk[Ag,Ta]w^ - AgTo^aW + AgTd^yjO + AgRjdkW, a)^ - dkT^^^a . 

R\ Rl Rl ^ S| ' 

By using the fact that: 

9+4 

Rl= Yl dk[Aq,S^>^^a]A^>w, 

q'=q-i 

Using the definition of the operator Aq leads to: 

[Aq, S^>_^a]A^>w{x) = - j h{y){S^,_^a{x) - S^,_^a{x - 2-'iy))A^,w{x - 2"^y)dy. 
and: 

\[Aq,S^>_^a]A^,w{x)\ < \\VS^>_^a\\Loo2-'^ j 2'^^ \h{2'iu)\\2'iu\\A^,w\ix - u)du, 
<2'^''\\VS^,_,a\\Lo.\{h{2'i-)\-\*A^,w)\{x). 

So we get: 

\\[Aq,S /_^a]A /w\\lp < \\VS '_-^^a\\L°o\\A 'w\\lp 
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we readily get under the hypothesis that a < 1, 

'^""WRlhr' < Yl '^"'\\^Sg>_^a\\L^\\Ag>w\\LP. (8.58) 

9 =9-4 

We have then: 

2'"^P^IIl. < c^llVall^.-^ 11^11^.+!-.. (8.59) 
In the case a = 1, we get: 

'^'^"WRIWlp <c„\\Va\\j.o \\w\\p.a+i-c (8.60) 

For bounding R^, standard continuity results for the paraproduct insure that if a < 1, 
Rq satisfies that: 

2^'^||i2^|kp < Co||Va||oa-l oo-+l-a. 

and if a = 1 

2^'^||i?^|kp < c„||Va||oa-i||tf||oa+i-a. 

II y II "^oo,l "-Op.r 

Standard continuity results for the paraproduct insure that R^ satisfies: 

2^"II-R'IU. <c,||V.i;|| .__^||a|| (8.61) 

-'-'00,00 ^p,r 

provided cr — a — y < 0. 
li(j-a-^ = then: 

2^"II^^IU-^c,||VT/;|bo^J|a|| (8.62) 

Next, standard continuity result for the remainder insure that under the hypothesis 
(7 > — -J, we have: 

'^'''\\Rt\W<c,\\Vw\\^.-4a\\ (8.63) 
For bounding R^ we use the decomposition: 



B ^ 



^1=11 5fe('5,'+2^9^\'«)> 
(?'>g-3 

which leads (after a suitable use of Bernstein and Holder inequalities) to: 

< ^ 2(^-«')("+f-i)2^('^+i-")||A,^|U.2«'(f+")||A^,a|U.. 



q>q-2 



Hence, since a + ^ — 1 > 0, we have: 

2'i''\\Rl\\Lv <c„||V 



LP >, Cq\\\/ W\\r,<T + l-a\\a\\ N_ 
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Combining this latter inequality with ()8.59l) , ()8.6ip and ()8.63p , and using the embedding 
Bp,r ^ -Boo.cx) for r = ^ + a — 1,(7 — a completes the proof of (18.560 . 

The proof of (18.570 is almost the same: for bounding R^, R^, R^ and Rg, it is just 
a matter of changing into sup^. We proceed similarly for R^. □ 
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